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Abstract 

It is usually accepted to consider an apsidal motion in binary stars 
as a direct confirmation that a substance inside stars is not uniformly 
distributed. It is shown in this paper that the apsidal motion in bi- 
nary systems observation data is in a good agreement with an exis- 
tence of uniform plasma cores inside stars if they consist of hydrogen- 
deuterium-helium mixture. 

PACS: 64.30. +i; 95.30.-k; 97.10.-q 



An apsidal motion (an advance of periastron) of close binary stars is a result 
of their non-keplerian moving. It originates from a non-spherical form of 
stars. The last has been produced by a rotation of stars or their mutual 
tidal effects. The traditional approach to the estimation of a value of this 
effect needs to suppose some redistribution of mass inside stars to reach 
an accordance between the measured data and calculations. Usually it is 
necessary to assume that the density of a substance at the central region of 
a star is hundred times greater than a mean density of a star [|l], [|2|. 

On the face of it, the model of constant parameter where den- 

sity inside a star core is constant, is in direct contradiction with the data 
of observation of an apsidal motion in close binary systems However, 
the traditional solution of the task about an apsidal motion is obtained for 
electrically non-charged matter. It is obvious that for electrically charged 
matter this solution must be reconsidered. 
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Let us estimate a contribution of inherent rotation of stars in an apsidal 
motion. According to 0-0, the ratio of u the angular velocity of an apsidal 
motion and Vt the angular velocity of rotation of a star about its axis is 

u ^?, {Ia- Ic) 

Vt 2 Ma2 ' ^ ' 

where I a and Ic are momenta of inertia relative to principal axes of the 
ellipsoid: 

lA-Ic = ^{a'-c% (2) 

where a and c are polar and equatorial radii of stars. 
Thus 

a; _ 3 (g^ - c^) 

" 10 a? ■ ^'^> 
It is correct when an angular velocity on the ellipse coincides with VL an 
angular velocity of rotation of a star about its axis . 

In the absence of rotation, the equilibrium equation of plasma inside a 
star is [H: 

IQg + PgEg = (4) 

where qg, Pg and Eg are acceleration of gravitation, gravity-induced 
density of charge and intensity of gravity-induced electric field, {div qg = 
47r G 7, div Eg = '^t^Pg and pG = y/G'y). 
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One can suppose that at rotation, an additional electric charge with density 
Pq and electric field Eq can arise under action of a rotational acceleration 
gn, and the equihbrium equation takes the form: 

(7g + 7n){gG + gn) = {pG + pn){EG + En), (5) 

where 

div {Eg + En) = 47r(pG + Pn) (6) 
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or 



div En = Airpn. (7) 

We can seek a solution for electric potential as an expansion on spherical 
functions. Assuming that an eccentricity is small, we can limit by a second 
term of the expansion: 

ip = Cn r\3cos'^e - 1) (8) 
or in Cartesian coordinates 

^ = Cn(3^'-x'-y'-^'), (9) 

where Cq is a constant. 
Thus 

E,^2Cnx, Ey^2 Cn y, E, ^ -A Cn z (10) 

and 

div En = (11) 
we obtain important equations: 

Pn = 0; (12) 

^gn = pEn. (13) 
Since a centrifugal force must be counterbalanced by the electric force 

7 Q2 X = p 2Cn X (14) 

the constant is 

The potential of a positively uniformly charged ball is 



g/3 



K,2 



A negative surface charge induces inside the sphere an electric potential 



V{R) = -% (17) 



where according to [§[ Q = v GM, M is the mass of a star. 
Thus, the total potential inside the star under consideration is 



^. = ^ +^'-=(3c-^»-l) (18) 



2R [ ~W)^2ic 
The electric potential must be equal to zero on the surface of a star, since 
at = and at r = i? we obtain 

R'-Rl = ^,, (19) 

3 ^ 

where Ro is the radius of a star without rotation. This equation describes 
the form of a rotating star. From this equation 

a'-c' = (20) 

and Eq.(l) gives 

- = (21) 

n 407rG7 ^ ' 

When both stars of a close pair induce an apsidal motion, this equation 
transforms to 

^ = IL^(L + 1\ (22) 

n 407rGV7i 72/' ^ ^ 

where 71 and 72 are densities of stars. 
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A density of non-degenerate non-relativistic plasma can be determined in 
the following way. At high density and temperature an electron gas is al- 
most ideal |Q and the ideal gas law will serve as its equation of state in first 



4 



approximation. It needs to take to account two corrections to describe elec- 
tron gas properties more exactly - a correction for identity of electrons and 
a correction for a presence of positively charged nuclei. The first of them is 
positive, it increases the incompressibility of electron gas. The second cor- 
rection decreases its incompressibility and it is negative. These correction 
are known P],and that is why with allowance for both corrections the free 
energy of electron gas takes the form 

where Z is the nucleus charge, Oq = ^^1^ Bohr radius. 

At a constant full number of particles in the system and at constant 
temperature, the free energy minimum equation is 

what allows one to obtain the steady-state value of density of particles of 
hot non-relativistic plasma 

n, = ~ 2 ■ lO^^Z^ cm-\ (25) 

97r^ao 

Thus, the equilibrium density of a star is 
where A is the mass number, rrip is proton mass. 



If we introduce -P = the period of ellipsoidal rotation and = 77 the 
period of the apsidal motion, we obtain from Eq. (p2D the main equation of 
our theory: 



5 



where the constant 



160 G nip 

For hydrogen the value is equal to 1, for deuterium 0.5, for helium 

1/128. Thus, for two stars both consisting of hydrogen the right part of 
Eq.(]27|) is equal to 2. For pairs consisting of hydrogen-deuterium it is 1.5. 
For D-D pairs and H-He pairs it equals to 1. For D-He pairs 0.5, and for 
He- He pairs 1/64. 

Thus, for pairs of a star consisting of hydrogen, deuterium and helium 

The periods U and P are measured for few tens of stars and we can 
compare our calculation with the data of these measurements. The distribu- 
tion of binary stars by the value of P^^/U is shown on Fig.|I| in logarithmic 
scale. On upper abscissa the value of iog{^^ + -j^) is shown. The lines 
mark the values of parameters for different pairs of binary stars. 

It can be seen that all measured data (with the exception of few pairs only) 
are limited by factors 1 and 1 / 64. This argues for adequate interpretation and 
satisfactory accuracy of our estimations of the effect of the apsidal motion. 

The interesting result can be obtained if we additionally take in to account 
the measured data on mass of rotating pairs. Accordingly 0], the star mass 
is determined by the A/Z ratio: 



A _ 12Mq 
Z ~ V 

and from Eq.BS 



(30) 



AZ^ = (31) 
16mp ^ ^ 

At solution of this equations with using of data observation of and 
7^ from the apsidal motion we can obtain a distribution of stars on Z and 
A separately (Fig.^). As the periods of the apsidal motion of stars consist- 
ing of nuclei heavier than helium are more than 10^ years and rather are 
immeasurable, this distribution is cut off on 1 < Z < 2 and 1 < A < 3. 
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Figure 1: The distribution of binary stars on value of P^^/U. 
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Figure 2: The distribution of binary stars on values A and Z. 



This distribution rather does not give a new astrophysical information, but 
it illustrates an efficiency and a cardinal utility of developed consideration. 
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